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Abstract. Two spectral regimes of magnetic field amplification in magnetohydrody- 
namic (MHD) flows can be distinguished by the scale on which fields are amplified 
relative to the primary forcing scale of the turbulence. For field amplification at or 
below the forcing scale, the amplification can be called a "small scale dynamo." For 
amplification at and above the forcing scale the process can be called a "large scale 
dynamo." Non — local (in wave number) effects play a key role in both the growth 
of the small scale field in non-helical turbulence and the growth of large and small 
scale fields in helical turbulence. Mean field dynamo (MFD) theory represents a simple 
semi-analytic way to get a handle on large scale field amplification in MHD turbulence. 
Helicity has long been known to be important for large scale, flux generating, externally 
forced MFDs. The extent to which such MFDs operate "slow" or "fast" (dependent or 
independent on magnetic Reynolds number) has been controversial, but there has been 
recent progress. Simulations of dynamos in a periodic box dynamo and their quench- 
ing can now be largely understood within a simplified dynamical non-linear paradigm 
in which the MFD growth equation is supplemented by the total magnetic helicity 
evolution equation. For dynamos, the large scale field growth is directly related to 
the large scale magnetic helicity growth. Magnetic helicity conservation then implies 
that growth of the large scale magnetic helicity induces growth of small scale magnetic 
(and current) helicity of the opposite sign, which eventually suppresses the a effect 
driving the MFD growth. Although the o? MFD then becomes slow in the long time 
limit, substantial large scale field growth proceeds in a kinematic, "fast" phase before 
non-linear asymptotic quenching of the "slow" phase applies. Ultimately, the MFD 
emerges as a process that transfers magnetic helicity between small and large scales. 
How these concepts apply to more general dynamos with shear, and open boundary 
dynamos is a topic of ongoing research. Some unresolved issues are identified. Overall, 
the following summarizes the most recent progress in mean-field dynamo theory: 

For a closed turbulent flow, 

the non-Unear mean field dynamo, 

is first fast and kinematic, then slow and dynamic, 

and magnetic helicity transfer makes it so. 



1 Small Scale vs. Large Scale Field Amplification 

A dynamo is a process which exponentially amplifies or sustains magnetic energy 
in the presence of finite dissipation. In this paper I will focus on magnctohydro- 
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dynamic (MHD) dynamos, where the only flux dissipating term in Ohm's law 
for the total magnetic field is the resistive term. 

The simple definition of a dynamo given above does not distinguish the scale 
on which the magnetic energy is sustained against turbulent forcing, whether a 
net flux is produced in the spatial region of interest, or the nature of the forcing 
(e.g. shear driven or isotropically forced). It is helpful to distinguish between 
"small scale dynamos" which describe field generation at or below the turbulent 
forcing scale, and "large scale dynamos" which describe field generation on scales 
larger than the forcing scale. Both present their own set of problems. 

For the Galaxy, supernovae dominate the nearly isotropic turbulent forcing, 
and although there is a range of forcing scales, typically the dominant scale is 
~ 50 — lOOpc A cascade leads to a nearly Kolmogorov turbulent kinetic 
energy spectrum. Faraday rotation and synchrotron polarization observations 
reveal the presence of a random component of the Galactic field, also with a 
dominant scale of 50 — 100 parsecs, and an ordered toroidal field on the scale 
of > Ikpc 

An important point about large scale field growth is that regardless of long 
standing debates about whether the large scale fields of the Galaxy is primordial 
or produced in situ, [H,||J|] and regardless of similar debates about the origin of 
large scale, jet producing, poloidal fields in accretion disks (e.g. (5)) one thing 
we do know is that the Sun, at least must have a large scale dynamo operating 
because the mean flux reverses sign every 11 years. If the flux were simply flux 
frozen into the sun's formation from the protostellar gas, we would not expect 
such reversals. 

For both small and large scale dynamos, the richness and the complication of 
dynamo theory is the non-linearity of the MHD equations. To really understand 
the theory, we need to understand the backreaction of the growing magnetic 
fleld on the turbulence, the saturated spectra of the magnetic field, and the 
spectral evolution time scales. Do theoretical and numerical calculations make 
predictions which are consistent with what is observed in astrophysical systems 
or not? What are the limitations of these predictions? 

1.1 Small Scale Dynamo 

The "small scale dynamo" describes magnetic field amplification on and below 
the turbulent forcing scale (e.g. I^j^j^). Field energy first builds up to equipar- 
tition with the kinetic energy on the smallest scales because the growth time 
is the turnover time, the turnover times are shorter for smaller scales, and the 
equipartition level is lower for the small scales. The fast growth, seen in non- 
linear simulations, can be predicted analytically |^,0,||,D. The approach to near 
equipartition is not controversial but the shape of the saturated small scale spec- 
trum needs more discussion, particularly for the Galaxy. 

Recent simulations of forced non-helical turbulence for magnetic Prandtl 
number (Pr = Vy/vM, where Vy is the viscosity and vj^i is the magnetic dif- 
fusivity) satisfying Pr > 1 in a periodic box have shown that the field does not 
build up to anywhere near equipartition on the input scale of the turbulence 
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[ P^IPH - Rather, the field piles up on the smallest scales. The interpretation is 
that the forcing scale velocity directly shears the field into folds or filaments with 
length of order the forcing scale but with cross field scale of order the resistive 
scale, accounting for the spectral power at small scales, (see for a careful 
semi-analytic study of the Pr >> 1 case). Since the small scale field is amplified 
by shear directly from the input scale, one can argue that there is a non — local 
direct cascade. 

How does the observed small scale field of the Galaxy compare with the 
above results? Since observations show near equipartition of the field energy 
with the kinetic energy at the forcing scale there is a discrepancy between 
the observations and the numerical simulation results for non-helical turbulence. 

A possible resolution might arise from the idealized problem, explored in 
[ p!4| . When the turbulence is forced with sufficient kinetic helicity (v • V x v), 
the spectrum changes in an important way. Figures 1 and 2 below come from 
3-D MHD simulations in which kinetic energy is forced in a box, and the helical 
fraction of the kinetic forcing is varied {fh = 1 corresponds to maximal helical 
forcing). The forcing wave number is k/2Tr = 4.5. For sufficiently large fh, (de- 
termined by that for which the kinematic dynamo can grow) the magnetic 
spectrum grows two peaks, one at A:/27r = 1 and one at the forcing scale. The 
kinetic helicity thus influences both the small and large scale field growth. 

Whether these results apply to the Galaxy or protogalaxy is unclear since 
those systems have shear (see also jl^), boundaries, and stratification unlike 
the simulations of . But the principle that the large and small scale fields are 
both infiuenced by helicity has been demonstrated. 

The emergence of two peaks motivates a two-scale approach which is a great 
simplification that will be exploited later. The shift of the small scale peak to 
the forcing scale is less understood than the generation of the large scale field. I 
now discuss the latter. 

1.2 Large Scale Dynamo 

The large scale field of the Galaxy seems to be of quadrapole mode and 

therefore the planar component of the large scale field has the same sign across 
the mid-plane. The toroidal field reverses on scales of a few kpc in radius In 
these annuli, there appears to be a net toroidal flux when integrated over the full 
height of the Galaxy. If this inference continues to survive future observations, 
the mechanism for field production must produce a net toroidal flux in annuli 
that extend the full vertical disk thickness, not just a net magnetic energy. 

A leading framework for understanding the in situ origin of large scale mag- 
netic field energy and flux growth in galaxies and stars, and even for the peak 
in the large scale magnetic energy in the helically forced case of Figure 1 (see 
also |l^) has been the mean field dynamo (MFD) theory M,|8|,p7|Jl8| . The the- 
ory appeals to some combination of helical turbulence (leading to the a effect), 
differential rotation (the f2 effect), and turbulent diffusion (the /? effect) to ex- 
ponentiate an initial seed mean magnetic field. Ref. [p^ developed a formalism 
for describing the concept |Q that helical turbulence can twist toroidal ((/)) fields 
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Fig. 1. Saturated kinetic and magnetic energy spectra for successive values of frac- 
tional helicity fh- 
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Fig. 2. Time sequence of kinetic and magnetic energy spectra for fh = l 
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into the poloidal (r, z) direction, where they can be acted upon by differential 
rotation to regenerate a powerful large scale toroidal magnetic field. The turbu- 
lent diffusion serves to redistribute the flux so that inside the bounded volume 
of interest, a net flux can grow. 

The formalism separates the total magnetic field into a mean component B 
and a fluctuating component b, and similarly for the velocity field V. The mean 
can be a spatial or ensemble average. The ensemble average is approximately 
equal to the spatial average when there is a scale separation between the mean 
and fluctuating scales. In reality, the scale separation is often weaker than the 
dynamo theorist desires, though a weak separation is also helpful given the 
limited dynamic range of simulations. I proceed to consider spatial averages 
to simplify the discussion. (Please also note that non-helical large scale field 
generation exists from the magnetorotational instability (MRI) but I do not 
consider this here and focus on externally forced systems. Actually the helical 
MFD dynamo may be operating even in systems with the MRI. Also I do not 
consider the model of here as that will be covered elsewhere in this volume.) 

The mean field B satisfies the induction equation [T^ , p^ 

- = -cVxE, (1) 

where 

E = - (V X B) /c - (v X b)/c + i/M V X B, (2) 

(v X b)j = a.jBj - l3^jkdjBk (3) 

is the turbulent electromotive force (EMF), and vm — rjc^/AiT is the magnetic 
diffusivity defined with the resistivity rj. Here aij contains Parker's twisting 
(the a effect) and Pijk contains the turbulent diffusivity. Ref. calculated 
E to first order in B for isotropic Uij and Pijk and hence the pseudo-scalar 
and scalar dynamo coefficients a and (3 respectively to zeroth order in B from 
the statistics of the turbulence, ignoring the Navier-Stokes equation. When the 
Navier-Stokes equation is not used, we speak of the "kinematic theory." Using 
the equation for the fluctuating field dth, plugging it into (v x b), the standard 

approach gives a kinematic a oc -J{v{t) ■ V X v(t'))dt' Tc(v • V X w)dt and 

/3 oc / (v(i) • v{t'))dt' ~ Tc(v • v)dt. I will come back to correcting the form for 
a in section 2.3 and 2.4. because we really want a fully dynamic theory, that 
accounts for the dynamo coefficients' dependence on b and B. Only then can 
one fully address the fundamental problem of mean field dynamo theory: how 
does the growing magnetic field affect the rate and saturation level of growth? 
Substituting (||) into (|l|), gives the MFD equation 

_==Vx (VxB)+Vx (aB)-Vx (/3-|-A)VxB. (4) 
at ' 

The first term on the right is the non-controversial "J7 — effect." If one assumes 
V, a and /3 to be independent of B, then (||) can be solved as a linear eigenvalue 
problem for the growing modes of B in the Sun and other bodies. However 
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a rapid growth of the fluctuating field necessarily accompanies the MFD in a 
turbulent medium. Its impact upon the growth of the mean field, and the impact 
of the mean field itself on its own growth have been controversial. 

The controversy arises because Lorentz forces from the growing magnetic field 
react back on the turbulent motions driving the field growth [p| j2l] , |2^ , p3| , ^4|j25|p7[ | . 
It is useful to distinguish between fast MFD action (also called "rapid" MFD 
action) and, slow MFD action (also called "resistively limited" MFD action). 
Fast MFD action proceeds at growth rates which do not go to zero as the mag- 
netic Reynolds number Rm —> oo, and maintain this property even when the 
non-linear backreaction from the magnetic field is included. Slow MFD action 
proceeds at rates that vanish as Ri\i — > oo. I sometimes use "resistively limited" 
rather than "slow" because the former more explicitly describes the reduced 
action. 

For galaxies and stars, conventional wisdom (which could be challenged |p8| 
and see section 2.3) presumes that rapid MFD action is necessary if the observed 
large scale fields are to be wholly produced and sustained by the MFD, given 
observed cycle periods and available time scales. That this may not be the case is 
a separate issue from understanding what the theory can actually provide. The 
latter is the focus herein. If dynamos in stars and galaxies do operate fast, then 
we would like to understand how, in light of recent numerical and theoretical 
evidence for slow dynamos in periodic boxes . What are the differences 

between these dynamos and real systems? where does the theory fail and where 
does the theory succeed? 

There are several issues which must be disentangled. First is the role of 
magnetic helicity conservation in constraining dynamo theory. In the steady 
state, such constraints are strongly infiuenced by boundary conditions, so one 
must be careful to understand the differences when applying idealized equations 
for closed systems to real systems with boundaries. Time-dependent dynamical 
constraints require helicity conservation to be supplemented in some way by 
the Navier-Stokes equation. Then the fully dynamical evolution of a and the 
non-linear backreaction can be studied. 

Two directions emerge. One is to produce a time dependent dynamical non- 
linear theory that fully agrees with numerical simulations in periodic boxes. 
This has recently been done The second is to recognize that in real sys- 
tems, there is physics which has not yet been comprehensively studied. This 
includes shear |l5y29|] , boundary terms which alleviate magnetic helicity conser- 
vation constraints, gravity, and a vertical variation of the dynamo coefficients. It 
has proven difficult in the simplest generalizations, to make the dynamo asymp- 
totically fast |^,|9j, but this may not be needed in some applications (e.g. 
Galaxy) if a kinematic phase lasts long enough. This will be discussed in detail. 

In section 2 I discuss the role of magnetic helicity conservation in dynamo 
theory. I first show that the escape of magnetic helicity through the boundaries 
might play an important role in maintaining fast MFD action during the steady 
sustenance phase of dynamos, and thus in interpreting quenching studies for this 
regime. I discuss the direct observational implications of the boundary terms and 
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the fact that stars and disks harbor active coronae. I then show that the time 
dependent mean field dynamo is really a process by which magnetic helicity 
gets transferred from small to large scales by a non-local inverse cascade and 
that the time-dependent dynamical quenching in recent periodic box simulations 
can be understood in this framework At late times, these box dynamos 
are resistively limited, depending sensitively on the magnetic Reynolds number, 
however the kinematic phase lasts quite awhile and this fact has very important 
implications. In the last part of section 2, 1 discuss an unsolved puzzle that arises 
in the derivation of the successful dynamical quenching model. In section 3, I 
summarize some key conclusions and pose open questions. 

2 Magnetic Helicity Conservation and Dynamo 
Quenching 

Although the MFD theory predates detailed studies of magnetohydrodynamic 
(MHD) turbulence, the MFD may be viewed as a framework for studying the 
inverse cascade of magnetic helicity. Whether this inverse cascade is primarily 
local (proceeding by interactions of eddies/ waves of nearby wave numbers) or 
non-local (proceeding with a direct conversion of power from large to small wave 
numbers) is important to understand. The simple MFD seems to be consistent 
with the latter [l6| . 

From the numerical solution of approximate equations describing the spectra 
of energy and helicity in MHD turbulence, Ref. ||3l| showed that the a effect 
conserves magnetic helicity, Hm ~ /(A-B)d'^x, by pumping a positive (negative) 
amount to scales > / (the outer scale of the turbulence) while pumping a negative 
(positive) amount to scales ^ /. Magnetic energy at the large scale was identified 
with the B of Ref. Thus, dynamo action leading to an ever larger B, hence 
the creation of ever more large scale helicity, can proceed as long as helicity of the 
opposite sign can be removed or dissipated. More on this in section 2.3. Here I 
first derive the general magnetic helicity conservation equations used to constrain 
the turbulent EMF, and investigate the implications of two steady state cases in 
detail before considering the time dependent case and interpretation of dynamo 
simulations. 

Using Ohm's law for the electric field, 

E = -c^^V X B + 77J (5) 

and averaging its dot product with B, gives 

(E - B) = E - B + (e-b) = -c'^v x b) • B + TyJ • B + (e • b) (6) 

where J is the current density. A second expression for (E-B) also follows from 
Ohm's law without first splitting into mean and fluctuating components, that is 

(E • B) = (J • B) = 7;J • B + 7y(j ■h)=i]l-B + c'^vm (b • V x b) . (7) 
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Using and (^), we have 

-c-\v xh) ■B = c-^X (b • V X b) - (e • b), (8) 

which can be used to constrain (v x b) in the mean field theory. 

Now consider E in terms of the vector and scalar potentials A and ^: 

E = -V<?- (l/c)atA. (9) 

Dotting with B = V x A we have 

E • B = - V<?> • B - (l/c)B • dtA. (10) 

After straightforward algebraic manipulation and application of Maxwell's equa- 
tions, this equation implies 

E • B = -(1/2)V • <PB + (1/2)V • (A X E) 

-{l/2c)dt{A ■ B) = (-l/2c)a^ff^(B) = ?7J • B, (11) 

where 



H^'{B) = {Ho, Hi) = [A • B, cm - cA x E] (12) 

is the magnetic helicity density 4- vector |Q , and the contraction has been done 
with the 4x4 matrix 77^^ where 77^^ — for fi ^ v, rj^^^ — 1 for fi = v = and 
Ty^^ = — 1 for /i = > 0. Taking the average of (^l]) gives 

d^H''{B) = -2c(E • B) = -2cE • B - 2c(e • b) = -2cr;(J ■ B). (13) 

If, instead of starting with the total E as in (||) , I start with e and then dot 
with b and average, the analagous derivation replaces (^ by 

5^ff''(b)--2c(e.b), (14) 

where H^ ih) indicates the average of _ff^(b). The latter is defined like ( p^ 
but with the corresponding fluctuating quantities replacing the total quantities. 
Similarly, starting with E and dotting with B, gives 

9^ff^(B) = -2cE-B, (15) 

where iJ^(B) is defined as in ( p^ but with the corresponding mean quantities 
replacing the total quantities. 

In the next three subsections, I consider the implications of magnetic helicity 
conservation for 3 separate cases. In the first I assume a steady state and ignore 
boundary terms. In the second, I include boundary terms but still demand a 
steady state. The third is the case in which boundaries are ignored, but a full 
time evolution of a is considered and a fully dynamical theory is obtained that 
agrees with recent numerical simulations. 
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2.1 Case 1: Homogeneous, Stationary Turbulence in Periodic Box 

Consider statistically stationary turbulence, where the averaging is over periodic 
boundaries. Then, the spatial divergence terms on the left of ( ]l3| ) become surface 
integrals and vanish. Discarding the spatial divergenece terms in (llj) gives 



(e.b) = -lai(a.b). 



(16) 



In the steady state, (|T^ vanishes and eqn (||) shows that (v x b) is resistively 
limited. 

If a uniform B is imposed over the periodic box, then B cannot change with 
time, and has no gradients. This is the case of Ref. which measures a but 
cannot test for MFD action. No mean quantity varies on long time scales. In this 
case, (ph then implies 



- c-i(v X b) -B = aB /c = c^^A (b • V x b) , 
where a — 033 for a uniform field in the z direction. Rearranging gives 

A(b- V X b) 



a 



B 



We then have 



l«l< 



k2b2T2 




ao 









(17) 



(18) 



(19) 



where 62 and V2 are the fluctuating field magnetic energy of the dominant en- 
ergy containing eddies (which is the forcing scale for both v and b for fh = 1 
turbulence as described in section 1), k2 is the wavenumber for that scale (= 
the forcing scale), and n = 3/4 if the current helicity is dominated by large 
wavenumber and n = 1 if it is dominated by small wavenumbcrs. The quantity 
T2 is the associated eddy turnover time and Rm.2 — 'Vilk2VM is the magnetic 
Reynolds number associated with the forcing scale. Assuming a steady state, I 
used V2 ~ &2 which is roughly consistent with numerical results. Then assum- 
ing forcing with maximal kinetic helicty, I replaced the numerator with ag, the 
maximum possible value of a helical quantity of that dimension. 



Pouquet Correction and Connection to Previous Studies: If we now 

take a to be of the form first derived by js^, (discussed further in section 2.3) 
in the context of a maximally helical, force free o? dynamo in a periodic box, 
we have 



a = -(l/3)r((v • V X v) - (b ■ V X b)). 



(20) 
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where t is the correlation time of the turbulence at k2 ■ Consider the implications 
of this formula in the steady state. If we use (|l|) in (|), and allow for a non- 
uniform mean field it is straightforward to show that a can be written 



ao + i?M,2/3(B • V X B)/B, 



a = - ' (21) 

1 + Rm,2B IBI^ 

Note that brackets around the current helicity are present because we allow for 
the fact that the mean field can have a scale smaller than the overall system 
(or box) scale. For example, the growth of the mean field at wavenumber fc = 1 
in a periodic box can occur when the fc = field is zero. Indeed, the limit of 
this equation for a mean field of zero curl produces exactly the result obtained 
numerically by result for uniform B. However, as is clear from this formula, 
it only emerges when B is uniform and the system is in a steady state. 

For steady-state but non-uniform B in the simplest case of a maximally 
helical dynamo in periodic box with no shear, the field energy and large 
scale helicity sustenance depend on E ■ B. If we assume for example that (3 = 
Po =constant (unquenched, one extreme limit) and ignore mean velocities, and 
use (HI) 

^(E . B) = a{B') - /3(B . V X B) = fl^l±^^^!li|Z^JlZ^ ( S ' ) - (3,{B • V x B) 

= ^ (#) ^ (B- V X B), 

(22) 

where B^q ~ f2 in velocity units. Notice that choosing a constant f3 — (3o in 
the steady state, leads to an EMF that is resistively quenched, and is the same 
as that with an artificially imposed symmetric, resisitive quenching of a and /?. 
These forms of a and (3 arc however misleading in the sense that although their 
combination is consistent with helicity conservation in the EMF, the division 
of a and P in this way is not the division which was conistent with our initial 
assumption of /? = /3o and so their forms are mere artifacts [^6| . Actually, in the 
saturated state, the current helicities of the large and small scale must be equal 
and opposite (this follows from ( [l3| ) with no divergence terms). We 

then have (J • B) ex a — ao- This implies, for our assumed /3 = /3o, and for large 
Rm 2 that 

« = 5 ■ (23) 

This is the actual steady state form of a when (3 is unquenched. Note that there 
is no Rm.2 depednence in the separate forms of a or /3 even though E • B is 
resitively limited: in the steady state (E • B) must satisfy (E • B) = VM-f ■ bbB). 
Ultimately, one can solve for B to see this demonstrated. 

If we were to instead assume a cc P from the outset, one finds 

(24) 



l + i?Af,2(a/ao + (s')/B2^ 
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Later we will see a deeper implication of the appearance of the current helicity 
in©. 

2.2 Case 2: Inhomogeneous Turbulence, Finite Boundary Terms 
and Implications for Coronal Activity in Steady State 

Now consider a system (e.g. Galaxy or Sun) to have volume << universal volume. 
Integrating ( [III ) over all of space, [/, then gives 

/ E-Bd^x^ -(1/2) / V ■^B(fix + (1/2) / V • (A X E) (fx 
Ju Ju Ju 

~-{l/2c)dt I A-BSx = -(l/2c)9tJ{(B) ^ [ t]^ ■ Bd^x, (25) 
Ju Ju 

where the divergence integrals vanish when converted to surface terms at infinity. 
I have defined the global magnetic helicity 

J{(B) = [ A B d^x, (26) 
Ju 

where U allows for scales much larger than the mean field scales. It is easy to 
show that a parallel argument for the mean and fluctuating fields leads to 

5t5{(B) =dt [ A • B d^x = -2c / E-Bd^x (27) 
Ju Ju 

and 

dM{h) =dt [ (a • b) d^x = -2c / (e . b) d^x = -2c / e-hd^x^ 5t J{(b), 
Ju Ju Ju 

(28) 

where the penultimate equality in (28) follows from redundancy of averages. 

I now split (|2^) and ( p8| ) into contributions from inside and outside the 
rotator. One must exercise caution in doing so because is gauge invariant and 
physically meaningful only if the volume U over which is integrated is bounded 
by a magnetic surface (i.e. normal component of B vanishes at the surface), 
whereas the surface separating the outside from the inside of the rotator is not 
magnetic in general. Ref. shows how to construct a revised gauge invariant 
quantity called the relative magnetic helicity, 

J{fl,,(B,) = J{(B„Po)-J{(P,;,P„) (29) 

where the two arguments represent inside and outside the body respectively, and 
P indicates a potential field. The relative helicity of the inside region is thus the 
difference between the actual helicity and the helicity associated with a potential 
field inside that boundary. The use of P^ is not arbitrary in ( p9| ) , and is in fact 
the field configuration of lowest energy. While (|2^) is insensitive to the choice 
of external field [M , it is most convenient to take it to be a potential field as is 
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done in ( p9D symbolized by Pq. The relative helicity of the outer region, JC/j.o, 
is of the form ( p9| ) but with the o's and i's reversed. The Jf^, is invariant even if 
the boundary is not a magnetic surface. 

The total global helicity, in a magnetically bounded volume divided into the 
sum of internal and external regions, U = Ut + Ue, satisfies 

J{(B) =Mfl,o(B)+J{i^,,(B), (30) 

when the boundary surfaces are planar or spherical. This latter statement results 
from the vanishing of an additional term associated with potential fields that 



would otherwise appear in ( |30[ ). Similar equations apply for B and b, so (^ 
and (ES) can be written 



5t?{(B) = ftJ{fl,o(B) + dt^RAB), (31) 

and 

dtMh) ^ dt^H^M + dt^R^h) (32) 
respectively. According to equation (62) of Ref. 25, 

dt^KR^B) = -2c [ E,-Bd^x + 2c [ (Ap x E) • dS, (33) 

JUi J DUi 

where Ap is the vector potential corresponding to a potential field V m.Ue, and 
DUi indicates surface integration. Similarly, 

^t'KR,^(^) = -2c / E ■ BSx + 2c [ (Ap x^ ■ dS (34) 



and 



dt'KR,i{h) ^ -2c I e hd^x + 2c ap x e ■ dS. (35) 

JUi JoUi 

Note again that the above internal relative helicity time derivatives are both 
gauge invariant and independent of the field assumed in the external region. If 
we were considering the relative helicity of the external region, that would be 
independent of the actual field in the internal region. 

Now if I take the average over a region < CA, I then replace (|3^) by 

(e • b) = -la,5{^_,(b) + (V • (ap x e)), (36) 

where the brackets indicate integrating over Ui or smaller. We now see that 
even if the first term on the right of (|3^) vanishes, (e • b) contributes a surface 
term to (||) that need not vanish. The turbulent EMF is not explicitly resistively 
limited as in the previous section, since the surface term can dominate. Thus in 
a steady state for Rm » 1, an outflow of magnetic helicity may enable fast 
MFD action. Since the right of ( [T^ ) is small for large Rm, the magnetic helicity 
flux has contributions from the small and large scale field. 
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Note that the sign of the pseudo-scalar a coefficient changes across the mid- 
plane of an astrophysical rotator. This means that the magnetic helicity gener- 
ation in one hemisphere is of opposite sign to that in the opposing hemisphere. 
One can then imagine that helicity flow could take place across the mid-plane; 
the loss of helicity from say the top hemisphere into the bottom, acts in congru- 
ence with the loss of the opposite sign of helicity from the bottom hemisphere 
into the top |^ allowing rapid MFD action. Even if there is not strong cou- 
pling at the mid-plane between the two hemispheres, an enhanced diffusion at 
the mid-plane side boundary of each turbulent/convective region could in prin- 
ciple dissipate the helicity require to allow a rapid MFD and a rapid net flux 
generation. 

But there are reasons why a fast MFD in stars and galaxies would more likely 
involve escape of magnetic helicity out of the external boundaries in addition 
to merely a redistribution across the mid-plane. First, buoyancy is important 
for both disks and stars and the Coriolis force naturally produces helical field 
structures in rising loops. Second note that the mean Galactic field appears 
to be quadrapole. This already requires some diffusion at the upper and lower 
boundaries of the disk 0. Third, the solar cycle involves a polarity reverse of the 
mean dipole field and this requires the escape of magnetic fields out through the 
solar surface. Fourth, note that for the Sun, the most successful dynamo models 
seem to be interface models where diffusion just below the base of the 
convection zone is reduced, with no alternative transport mechanism downward. 
Finally, note that some aspects of solar magnetic helicity loss have been studied 



and observed [B5 36 



Let us explore some implications of depositing relative magnetic helicity to 
a astrophysical corona. I assume that the rotator is in a steady state over the 
time scale of interest, so the left sides of (Q) and (^) vanish. For the sun, where 
the mean field flips sign every ~ 11 years, the steady state is relevant for time 
scales less than this period, but greater than the eddy turnover time (^^ 5 x 10^ 
sec). Beyond the ~ 11 year times scales, the mean large and small scale relative 
helicity contributions need not separately be steady and the left hand sides need 
not vanish. 

The helicity supply rate, represented by the volume integrals (second terms 



of Eqs. (34) and (|35|)), are then equal to the integrated flux of relative magnetic 
helicity through the surface of the rotator. Moreover, from ([T^), we see that the 
integrated flux of the large scale relative helicity, = 5'fl^i(B), and the integrated 
flux of small scale relative helicity, = ?'7?_i(b), are equal and opposite. We thus 
have 

'jR,i(B) = -2c I E ■ B(fix = 2c I e • hd^x = -Jfl„i(b). (37) 
Jui Jui 

To evaluate this, I use (0) and (|) to find 

E = -c"\aB-/3V X B), (38) 
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throughout Ui. Thus 

S'rAB) = -Jfl.,(b) = 2 / (aB' - /5B . V X B)d3x. (39) 

J (7, 

This shows the relation between the equal and opposite large and small scale 
relative helicity deposition rates and the dynamo coefficients. 

Now the realizability of a helical magnetic field requires its turbulent energy 
spectrum, E^^ , to satisfy |^ 

(b) > ^mkih)\ , (40) 

where !Kfc is the magnetic helicity at wavenumber k. The same argument also 
applies to the mean field energy spectrum, so that 

i?f (B) > ^fc|Jf,(B)|. (41) 

If I assume that the time and spatial dependences are separable in both E'^' and 
IK, then an estimate of power delivered to the corona can be derived. I presume 
that the change in energy associated with the helicity flow represents an outward 
deposition rather than an inward deposition. This needs to be specifically calcu- 
lated for a given rotator and dynamo, but since the source of magnetic energy 
is the rotator, and buoyancy and reconnection interplay to transport magnetic 
energy outward, the assumption is reasonably motivated. 
For the contribution from the small scale field, we have 

E'\h)^jE^Hh)dk>^Jk\Jk,RA^)\dk>^J\9kMAWk 



where the last equality follows from the first equation in (39). The last quantity is 
exactly the lower limit on i?*^(B). Thus the sum of the lower limits on the total 
power delivered from large and small scales is ^^|Ji^_i(B)| + |-F_R,i(b)| = 
2^^j^|If(B)|. Now for a mode to fit in the rotator, k > kmin — 27r//i, where /i is a 
characteristic scale height of the turbulent layer. Using (^9|), the total estimated 
energy delivered to the corona (=thc sum of the equal small and large scale 
contributions) is then 

> 2:^|:r;,,,(b)| = 2^|:r«,(B)| = X _ • v x b , (42) 

OTT oTT 11 ave 

where V is the volume of the turbulent rotator and ave indicates volume aver- 



aged. I will assume that the two terms on the right of (42) do not cancel, and 
use the first term of ( p2[ ) as representative. 

Working in this allowed time range for the Sun (5 x lO^sec < t < llyr) I apply 
( p2[ ) to each Solar hemisphere. Using the first term as an order of magnitude 
estimate gives 



\ 3 J V7xl0i"cmy V40cm/s; Vi50G; s' 
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I have taken a ~ 40 cm s^^ (a low value of a) and have presumed a field of 150G 
at a depth of lO^'km beneath the solar surface in the convection zone, which is 
in energy equipartition with turbulent kinetic motions Q]. 

As this energy is available for reconnection, Alfven waves, winds, and par- 
ticle energization, we must compare this limit with the total of downward heat 
conduction loss, radiative loss, and solar wind energy flux in coronal holes, which 



cover ~ 1/2 the area of the Sun. According to |40 , this amounts to an approx- 
imately steady activity of 2.5 x lO^^erg s~^, about 3 times the predicted value 
of (^3|). There is other evidence for deposition of magnetic energy and magnetic 
+ current helicity in the Sun |||,|6@ . 

Active galactic nuclei (AGN) and the Galactic interstellar medium (ISM) 
represent other likely sites of mean field dynamos j4^. For the Galaxy, E^'^ > 

(7ri?2)a# ~ 10''°(i?/12kpc)2 x{a/10^cm/s)(B /5 x 10-^G)2erg/s in each hemi- 
sphere. This is consistent with coronal energy input rates required by and 
For AGN accretion disks, the deposition rate seems to be consistent with 
what is required from X-ray observations. Independent of the above, the most 
successful paradigm for X-ray luminosity in AGN is coronal dissipation of mag- 
netic energy [Q. 

So boundary terms can potentially alleviate any helicity constraint and allow 
fast steady dynamo action. But how this specifically happens and the specific 
boundary physics of a real system will need more study to see what field wave 
numbers, if any, are preferentially shed. 



2.3 Case 3: Time-Dependent Dynamo Action and Dynamical 
Quenching in a Periodic Box 

For the question of actual field amplification and time scales ultimately associ- 
ated with cycle periods for a real dynamo, a time dependent, non-linear dynam- 
ical theory is required. Here the time derivative term of (Jl^) important. (Note 
also that time dependent dynamo like effects are also important for magnetic 
field adjustment in Reverse Field Pinches ji^ .) 

An important step forward was the work of Ref. based on the Eddy 
Damped Quasi-Normal Markovian spectral closure scheme. There an approxi- 
mate set of equations describing the evolution of magnetic and kinetic energy 
and helicity was derived. Refs. |lj,|l^ performed numerical simulations of the 
process, and Ref. further simplified the equations of by considering a 
two-scale approach for fully helically forced turbulence. In [|jO[ it was assumed 
that the large scale field grows primarily on scale fci and the small scale turbulent 
field is peaked at fc2 >> fci. The analytic model therein is largely consistent with 
the forced helical turbulence, periodic box simulations of Refs. [ p^Hl^ . It should 
be mentioned however, that while the large scale field dynamics are consistent 
with a single scale fci at all times, the wavenumber of the small scale peak seems 
to migrate from the resistive scale to the forcing scale, in which case the two-scale 
approach applies but with a migrating fc2 unlike that considered in jSOt . 

More work is needed to understand the migration of the small scale peak. For 
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present purposes we ignore this complication which does not effect the accuracy 
of the resulting fits of the large scale field growth. 

The basic concept of the successful dynamical quenching model Q for the 
o? dynamo is that the growth of the large scale field is the result of a segregation 
of magnetic helicity. Magnetic helicity of one sign grows on the large scale, while 
the opposite sign grows on the small scale, up to resistively limited conserva- 
tion. The growth of the small scale magnetic helicity also grows current helicity 
which suppresses a. The two essential equations needed in this analysis are (1) 
the equation for large scale magnetic helicity evolution (|l5|) and (2) the equation 
for total magnetic helicity evolution ([l3|). For more general dynamos with shear, 
the equations of the dynamical theory are (1) the vector equation for the large 
scale magnetic field, and (2) the equation for total magnetic helicity conserva- 
tion. In general, the paradigm that emerges is that the total magnetic helicity 
conservation acts as a supplementary dynamical equation that is coupled to the 
evolution of the large scale field equation. 

Ignoring the boundary terms in equation ([l^). and subscript 1 and 2 to 
indicate mean and fluctuating scales, and M to indicate "magnetic" we then 
define H^^ = iJ"(B)/2. Thus (|l|) implies 

dtH^' = 2a£;f - 2(A + /3) J • B, (44) 

where where = B^/2. Now we replace spatial derivatives of the mean scale 
with ki, and then note that the current helicity term on the right is related to 
the magnetic helicity J • B = /cJiJ*^. We then have 

dtHf = 2aEf - 2(A + l3)klH^^ . (45) 

Ref. describes the conditions for which this is consistent with the analogous 
equation that arises from the spectral treatment o f |pl| . 

The second equation we need is a re- write of (|l3[). Since the left- most side 
is the sum of contribution from large and small scales, as is the right- most side, 
we have in the two-scale approximation 

dtH^' + dtH^' ^ -2XklH^' - 2\klH^'. (46) 

We also need a prescription for a and (3. A configuration space approach for 
a and its pitfalls are described in the next section, but for now, let us extract 
the dynamical form derived from For a this is 

a = -(1/3)t((v V X v) - (b- V X b)), (47) 

where r is a correlation time of the turbulence at fc2- In externally forced sim- 
ulations, the first term on the right of ( |4^ ) (the kinetic helicity) is typically 
maintained at a fixed value. At early times, the second term on the right is 
small. Thus the second term on the right (current helicity) can be thought of 
as a backrcaction on the first term. It arises from inclusion of the Navier-Stokes 
equation. The idea that the second term represents a backrcaction was inves- 
tigated by for a stationary system and a correction to a was derived. 
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However, here I will derive the fully dynamical correction following ||3C 
re-write (E^ as 



"0 



Tfc2i?2 



M 



/3. 



First 



(48) 



If we assume that the kinetic helicity is forced maximally, and take r = 2/fc2f2, 
then ao = 2u2/3. 

Unfortunately, there is not yet convergence on a rigorous prescription for 
j3 in 3-D, but I will consider the two cases discussed earlier, /? = /3oa/ao a-nd 
P = Po = W2/fc2=constant as examples. In [p6| other prescriptions are considered, 
including those for which /3 is quenched but not resistively. In general, the need 
for the appropriate form of f3 is an important input, and determines ultimately 
the form of saturated a in the dynamo because the two are constrained 
through E • B as discussed at the end of section 2.1. However, the particular 
form of (3 is also less important for illustrating the role of helicity conservation 
in the dynamo than in dynamos with shear as for the former since a range 
of choices are all relatively successful. I will show that the difference for large 
^2/^1 between choices of (3 is really quite minimal. More prescriptions for (3 are 
considered in pGf . 

The same type of formalism which leads to the prescription for a, produces a 
prescription for f3 involving an approximate sum of kinetic + magnetic energies 
times a correlation time, at least to first order sort of motivating the case (3 = j3o. 
More work is ongoing and indeed for very strong mean fields since (3 must at 
least respond to the mean Lorentz force. The case (3 — (3oa/aQ is motivated by 
the fact that at late times one empirical combination of formulae which fit p6| 
has this relation. 

It is important to emphasize though, that as shown in section 2.1, a mislead- 
ing degeneracy emerges in the steady-state (or near steady state which amounts 
to the late time evolution) with respect to the choice of f3. Again, the reason is 
that for dynamos, it is really the (E • B) that matters for the growth of the 
large scale magnetic helicity, and thus the large scale field. The main effect of 
the prescription for /3 is the saturation value of the field. 

Following we need to solve (45) and (^). To do so, I rewrite them in 



dimensionless form. Define the dimensionless magnetic helicities hi = 2H^k2/v2 
and h2 = 2i7^^fc2/w| and write time in units of l/fc2i'2. I also define Rm = 
{v2/ki)/vM- (Note that this definition of Rm is based on the forcing-scale RMS 
velocity but on the large scale, k^^ . We will later employ the previously defined 
magnetic Reynolds number Rm,2 = RM{ki/k2)-) 

Using the above scalings we can replace (45) and ( ^6| ) with dimensionless 
equations given by 



and 



dthi = 



dth2 



4 fki 



3 Vfc 



hi{l + h2)-2hi 



ki 



k2RM 



kl 

k. 



2(1 



fc2 

/ll(l + /t2)+2/li-|(l 
"^2 



92^2) - 



92/12) 



2 h2k2 
Rm kl 



(49) 
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where (72 = in the above equations corresponds to P{t) — (3q ^constant, and 
52 = 1 corresponds to P{t) — a{t)(3o/ao. Solutions of these coupled equations 
are shown in Figs. 3-8, which are taken from |Q. The key parameters are ^2/^1, 
Rm, and (72- In the figures, I have compared these results to the empirical fits 
of numerical simulations in I used hi{t = 0) = 10"'^, but the sensitivity 
to hi(0) is only logarithmic (see ( [m] ) below). In Fig. 1, fc2/fci = 5 was used, 
following BOl, and in Fig. 2 fc2/fci =20 was used. 



In the figures, the solid lines represent numerical solutions to (49) and (|50|), 
whereas the dotted lines represent the formula given in ||l6|] , which is an empirical 
fit to simulation data assuming that a and /? are prescribed according to ( ^ 
and (|5^) below. More explicitly, Ref. |0 found that the growth of B was well 
described by the formula 

r2 / r2 



(1-B?/Sta*) 



. J-T r5 — 



where Bi q = Bi{t = 0). This can be rewritten using the notation above as a 
dimensionless equation for t in units of (^2^2)^^, namely 

h Lnpi/M(l-/iifc^/fc|)«^HfciA2~2/3)-i] 

2fci 2/3-fci/fc2 ■ 

Note that (|5^) and (|l]) correspond to a and /? quenching of the form 

" = i + s7b!M ^''^ 

and 

(3 = , , (53) 

where SB ~ i?A/(fci/fc2)(2/3-A;i/A:2) = i?M,2(2/3-fci/fc2), and i?M,2 = V2/k2VM- 
Eqns. (HI) and ( p3| ) are derived from those in by re-scaling Eq. (55) of 
with the notation herein. It can also be shown directly that, up to terms of order 
1/Rm, ( ^ ) is consistent with that derived by substituting ( p2|) and ( |53| ) into 
( p9| ) and solving for t. Note that in contrast to the suggestion of |16|, it is actually 
the forcing-scale magnetic Reynolds number, i?M,2, that plays a prominent role 
in these formulae. 

The solutions of (^) and ( |50| ) are interesting. Insight can be gained by their 
sum 

Q , , a , 2 f hiki ^ /l2fc2 



Rm V ^2 fci 



which corresponds to (46), the conservation of total magnetic helicity. If we 
make the astrophysically relevant assumption that Rm » 1, the right hand 



Recent Developments in Magnetic Dynamo Theory 19 



side of p^ ) is small for all hi and h2- It follows that dt{hi + /i2) =0 and for 
h(t = 0) = 0, this implies h2 = —hi. In this period, we can self-consistently 
ignore 1/Rm in (^)- If 92 = 1, this phase ends when /12 —1, so that /ii ~ 1. 
This is manifested in figure 5. 

This kinematic phase precedes the asymptotic saturation of the dynamo in- 
vestigated by other authors, in which all time derivatives vanish exactly. For 
this to happen, the right hand side of (|5^) must vanish, which is equivalent to 
demanding that /12 = —(^1/^2)^/11. Since the right hand sides of ( ^9|) and ( po|) 
are proportional to I + /12 when terms of order l/Rm are neglected, their vanish- 
ing requires that /12 = —1, and therefore, that hi = (fc2/fci)^. This is observed 
in figures 1 and 2. The asymptotic saturation (when the field growth ceases) 
takes a time of order tsat ~ RMk2/ki, which in astrophysics is often huge. Thus, 
although in principle it is correct that a is resistively limited (as seen from the 
solutions in figures 6 & 7) as suggested by |p3| , p5p^ , ^ , this is less important 
than the fact that for a time tkin < Rm the kinematic value of a applies. The 
time scale tkin here is given by a few kinematic growth time scales for the 
dynamo, more specifically, 

tMn - Ln[l//ii(0)](fc2/fci)/(4/3 - 2fci/fc2). (54) 

For hi{0) — 0.001, ^2/^1 — 5, tkin ^ 37, as seen in Fig 5. 

Note that tkin is sensitive to fc2/fci and independent of Rm- Figure 5 shows 
that there is significant disagreement in this regime with (^2]), but this formula 
was used in |jl^ only to model the regime t > Rm, so the result is not unexpected. 
We can see from the solution for a itself that indeed the solutions of ( p9| ) and 
(|o|) do match ( p^ ) for t > Rm (figures 6 and 7) . Figure 6 shows the difference in 
the a along with (^3|) for the two values Rm = 10^ and Rm = 10'^. Notice again 
the disagreement with the formula ( p2[ ) until t = Rm , and agreement afterward. 
This marks the time at which the resistive term on the right of (|4^) becomes 
competitive with the terms involving (1 -I- ^12). Asymptotic saturation does not 
occur until t ~ tsat = RMk2/ki as described above. 

Finally, note that q2 = corresponds to P — (3o- In general, this leads to a 
lower value of hi in the asymptotic saturation phase because this enforces zero 
saturation of /3, whereas there is still some saturation of a in this limit. (Note 
that q2 = corresponds to the case ofjB^ discussed in appendix B of Q.) For 
large fc2/fei the solutions of (^ ) and (|50|) are insensitive to (72 = or (72 = 1- This 
is because the larger fc2/fci, the smaller the influence of the q2 terms in (|9|) and 
(^0|). This is highlighted in figure 8 where the result for 52 = is plotted with 
the |l^ fit. This suggests that for large-scale separation, the magnetic energy 
saturation is insensitive to the form of P quenching. However, in real dynamos, 
magnetic flux and not just magnetic energy may be needed, so the insensitivity 
can be misleading because /3 is needed to remove flux of the opposite sign. From 
the low fc2/fci cases, it is clear that that (72 = 1 is a better fit to the simulations 
of@. 

The physical picture of the quenching process just described is this: helical 
turbulence is forced at fc2 (= 5 in p6[), and kept approximately constant by 
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Fig. 3. Solution for hi{t), fh = l,q2 = 1- Here k2/k\ = 5 and the three curves from left 
to right have Rm = 100, 250, 500 respectively. The dots are plotted from the formula 
used to quasi-empirically fit the simulations in Ref. 
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Fig. 4. Solution for h\{i), fh = 1, 92 = 1- Here ki/ki = 20 and the three curves 
from left to right have Rm = 10^, 10*, 10^ respectively. The dots are plotted from the 
formula for used to quasi-empirically fit simulations of Ref. mM 



forcing. Hence ao = —2tHY/^ ~ const. If Hi"^ , the magnetic helicity at ki 
(which reaches 1 here as a result of boundary conditions), is initially small — 
so that \2klH^ /S\ < |ao|, (|45|) (or (|49|)) shows that it wiU be exponentially 
amplified provided that the damping due to /3+ A does not overcome the a effect. 
Initially, a — ao, acting like a pump that moves magnetic helicity from /c2 to ki 
and driving the dynamo. This kinematic phase lasts until tkin as given by (^J). 
Eventually, the growing Hf^ results in a growing H^^ of opposite sign, which 
reduces a through H2 ■ i?Af-dependent quenching kicks in at t = tkin, but it is 
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Fig. 5. The early-time solution for hi{t), fh — 1, q2 = 1- Here for k2/ki = 5, and 
Rm — 10^, 10'', 10* from left to right respectively. Notice the significant departure 
from the formula of at these early times. For t < tkin there is no dependence on 
Rm and the growth proceeds kinematically. 



1 




Fig. 6. Solution of a/ao{t) for hi{t), fh — i, q2 = 0. Here fc2/fei = 5 and the solid lines 
are the solutions to and ( |50| ) for Rm ~ 10^ (top curve) and Rm = 10^ (bottom 
curve) respectively. The top and bottom dotted curves are from ([52[), interpreted from 
Ref. mA. Notice the long kinematic phase for the solutions, the overshoot, and the 
convergence with that of (M) at t = Rm for the Rm ~ 10^ case. 
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Fig. 7. This is the extension of the previous figure for later times. Notice the conver- 
gence of the Rm ~ 10^ solution to ( p2|) near t = Rm- 
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Fig. 8. Solution for hi{t), fh ~ 1, q2 ~ 0. Here k^/ki — 50 and the three curves from 
left to right have Rm = 10'^, 10^, 10* respectively. The dotted lines are plotted from 
the formula used to quasi-empirically fit simulations of Ref. |^^. For such large k2/k\ 
the fit to the data is only weakly sensitive to whether g2 = 1 or g2 = 0. 



not until t = Rm that the asymptotic formulae (^2|) and ( p3| ) are appropriate. 
Asymptotic saturation, defined by the time at which Bi approaches its maximum 
possible value of ik2/kjY/^V2, occurs at t = tgat = ^1/^2/^1- For t > Rm the 
numerical solution of (|49|) and (|50|), like the full numerical simulations of is 

53). The two-scale approach 
is also consistent with BOl in that magnetic helicity jumps from ^2 to fci without 
filling in the intermediate wave numbers. 

The emergence of the time scale tkm is interesting because it shows how 
one can misinterpret the implications of the asymptotic quenching formula ( p2|) 



well fit by the a in (p2[) with a corresponding /3 of 
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and ( |53| ) . These formulae are appropriate only for t > Rm ■ The large-scale field 
actually grows kinematically up to a value Bi = (fci/fc2)^/^W2 by t — tkin and 
ultimately up to Bi ~ (fc2/fci)^/^f2 hy t = tsat- For large Rm,2, these values 

1/2 

of Bi are both much larger than the quantity V2/Rj^ 21 which would have been 
inferred to be the saturation value if one assumed (B2f) and (m) were valid at all 



times. 

Dynamical quenching or time-dependent approaches recognizing the current 
helicity as a contributor to a have been discussed elsewhere [p|, ^9|j50| , ^ (see also 
[p5|), but here we have specifically linked the PFL a correction to the helicity 
conservation in a simple two-scale approach. Other quenching studies for closed 
systems such as [|5) and Q advocated values of a which are resistively limited 
and of a form in agreement with but with the assumption of a steady 



Bi. Assuming ([48D, and using (^5|) and (46) in the steady state, their formulae 
can be easily derived. However, one must also have a prescription for /3. If /3 is 
proportional to a, then the resistively limited formulae like (|24| ) emerges exactly, 
which is indeed approximately consistent with ( ^2|) and (|53| ) for large Rm- If 
j 3{t) = Po, as in then a formula for a without resistively limited quenching 
( P3|) emerges. On the other hand, fig. 8 shows that for large fc2/fci, the dynamo 
quenching is largely insensitive to (3. 

Interestingly, if we cavalierly apply these results for the Galaxy (by incor- 
rectly ignoring the shear and assuming an dynamo that produces force free 
large scale fields), and use k2/ki = 20, i?i(0) = 10~^, V2 = lOkm/s, and 
fc^^ ^ lOOpc, we would find the end of the kinematic regime to be at i = 140, 
or about 1.4 x lO^yr. After this stage the field growth would proceed very 
slow because of the large Rm, but the saturation field strength at this time 
is Bi ~ 1)2/4.5. Thus quite a large amplification can occur, even with an asymp- 
totically slow dynamo. However this is really an academic exercise since for the 
Galaxy we need to consider an a — ^2 dynamo, and the boundary terms. 

The generalized application of the principle that the magnetic helicity con- 
servations should supplement the mean field dynamo growth equation to account 
for the backreaction is applied more generally to dynamos with shear in [ p6[ . 

Finally, note that there is an important puzzle, hidden in the derivations 
here and those of [^^^ with regard to a quenching that we discuss in the next 
section. 



2.4 Deriving a in Configuration Space: a Puzzle 

The two-scale dynamical theory of section 2.3 for a-quenching is appealing 
because it nicely couples the equations and concepts of magnetic helicity evo- 
lution to the current helicity contribution in a, and fits simulation data well. 
The current helicity contribution was interpreted as a correction to the kinetic 
helicity contribution of kinematic theory. However, it depends on the current 
helicity contribution to a as presented in ( ^7|) being the total current helicity 
associated with ^2- To see what we mean by total and to show the complication, 
we consider the derivation of the turbulent EMF in configuration space. 
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The turbulent EMF can be written in three different ways: 
(v X b) = (v(0) X b(0)) + JoidMt') X hit'))dt' + /o*(v(t') X dtMt'M' 

- (v(t) X b(0)) + Jo (v(t) X dt'h{t'))dt'i55) 

= (v(0)xb(t))+/„*(9,,v(t')xb(t))di' 

The three hues in ( |55| ) simply correspond to the 3 relevant ways of using the 
formula f{t) = /(O) + df f{t')dt\ where g is an arbitrary function of time. 
If I assume that t >> 0, and that widely separated turbulent quantities do not 
correlate, the first terms on the right of the 2nd and 3rd lines respectively, can 
be dropped. We then have 

(v X b) = /J(v(t) X dt'h{t'))dt' = /o (5t'v(t') X h{t))dt'. (56) 



The second term in ( |56| ) as the standard textbook starting point |17 1^ for 



evaluating the turbulent EMF for a kinematic dynamo, but here I have not 
made any assumptions about the backreaction yet. 

To illustrate the point, consider the simple case in which VB — 0. The 
equation for the small scale field is then 

dt'h = B • Vv + b • Vv - V • Vb - V X (v X b) + i^A/V^b. (57) 



The penultimate term goes away when included in (|56| ) and we ignore the last 
term. The first terms on the right, upon the assumption that the dominant 
contributions to correlations are isotropic in v, gives the "textbook" expression 
for a plus extra terms, that is 



(v X b) ~ -5 f (v(t) . V X v{t'))dt' + Q{v^b). (58) 
3 Jo 

The terms symbolized by Q{v^b) are typically ignored using some version of the 
first order smoothing approximation. This is of questionable validity, given that 
the small scale field rapidly grows to exceed the mean field. We will come back 
to the relevance of these terms below. 

Now if instead I use the last term of ( p6| ) to expand the EMF, I must then 
invoke the Navier-Stokes equation for the time derivative of the turbulent veloc- 
ity 

dtv = -vVv-(v-Vv)-Vpe//+B-Vb+b-Vb-(b-Vb)+i/V2v+f(x,t), (59) 

where f is a forcing function and Pef f is the magnetic and thermal pressure. 
Upon plugging this into (|5^), the second and sixth terms on the right vanish. 
If we ignore the viscosity, and assume the dominant contribution to correlations 
are isotropic, we then have 

(v X b) ~ f / (b(i) • V X h{t'))dt' + Q{v\ fb) (60) 

3 .In 
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Notice two things about (|5^) and ( |60| ) : First they are equal to each other since 
they were derived from different choices of the expansion of (v x b). Second, they 
do not cleanly include the combination of the total residual helicity, required in 
when placed into (|5|) and (||) to derive (|9|) and @. The only way that 
(58) and ( |60| ) can have the form of the desired relative helicity is if Q{v^b) = 
f /o(b(t) • V X h{t'))dt' and if Q{v\ fb) = -f /o(v(t) • V x v{t'))dt'. I have 
been unable to prove that this is the case for the astrophysically relevant weak 
B regime. 

There is another approach to calculating the EMF in configuration space that 
does reveal a similar difference of helicities as that in (^), namely the approach 
of But V and b enter ( p7| ) whereas v^°) and b^°\ the statistically isotropic 
parts of V and b, enter |Q. To see this more explicitly, I write 

v = v(") +v('4), (61) 

where A indicates an anisotropic contribution, the result of the backreaction 
from B. Similarly, 

b = b(°)+b(^). (62) 

(Even when b is the result of stirring up an initial seed B, there is still a b^^-* 
which is the statistically isotropic part of b.) We then assume that the statistics 
of the zeroth order turbulent correlations are those of a homogeneous isotropic, 
"known" base state. The goal is to express turbulent correlations in terms of 
the zeroth order quantities. It is sufficient to demonstrate the basic idea invoked 
to all orders in B in with that derived to linear order in ||5^ Noting that 
(v X b)'-^-' vanishes, the lowest order contribution to the turbulent EMF is 

(v X b)(i) = (v(o) X b(i)) + (v(i) X b(0)) 
= loi^^^'Ht) X dt'h(^Ht'))dt' + /o(9t'v(i)(t') X b(o)(i))dt' ^ ' 

To linear order, using the induction equation for b^^^ and the Navier-Stokes 
equation for v'^^ it can be shown that by analogy to the derivations of ( p8| ) 
and (|60|), (combined with a revised first order smoothing approximation that 
assumes |b(-^^/B| < 1) ( |6^ ) becomes 

(v X b)(i) = • V X ^^°^{t'))dt' - ^*(b^"^(i) ■ V X b(")(i'))dt'^4) 

FBC showed that in the case of negligible mean field gradients, it is still the 
zeroth order kinetic and current helicities which appear most explicitly in a, 
even to all orders in B. 

It is clear that the zeroth order helicities are not necessarily equal to those 
constructed with the full turbulent quantities since 

(v V X v) = (v(") • V X v(")) + (vt-^) ■ V X v(°)) + (v(^) • V x v^"') (65) 
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and similarly 

(b • V X b) = (bW ■ V X b("'} + (b(^) • V X b(")} + (b^-^) • V X b("'}, (66) 

so the extra terms on the right must be dealt with. One might ask however, 
if the first terms on the right of ( |65|) and ( |6^ ) dominate, why can't we simply 
replace (b^^^ • V x b^^^) by (b • V x b) wherever the former occurs? The reason 
is that the appropriate helicity which then enters a is (b(°)-Vxb(o)) to lowest 
order. Then the H^^ entering on the left of ( ^5|) would be second order. But 
then the magnetic helicity, i?;^^, entering (^) would also be second order. Thus 
the current helicity entering a is zeroth order whereas that entering the helicity 
conservation equation would be second order. There is an ordering mismatch. 
This is a problem because the success of the model of 2.3 depends on our being 
able to circumvent this ordering ambiguity and presume that the current helicity 
entering p7| ) is exactly fc| times the ^ entering (^). 

The procedure outlined to derive (Q) and the subtlety just described with 
respect to ordering is basically the "ordering ambiguity" that was discussed in 
There it was shown that Refs. |^,^ effectively derived the form ( |63| ) 
rather than (^7|) by linearizing in terms of B but did not identify that they had 
derived the zeroth order contribution to a. Thus Refs. [ ^6p^ were actually using 
a similar expansion to that of The subsequent manipulations of p6| and 
required that they had derived a as a function of the full v and b, much 
like our manipulations in section 2.3. The issue also arises subtly in the k space 
derivation of PFL and is presently unresolved. 



3 Conclusions and Open Questions 
3.1 Small Scale Dynamo 

For non- helical turbulence, and for Pr > 1, current simulations indicate that the 
magnetic field piles up on the resistive scales k » kf when forced externally 
[ pT| in a periodic box. The reason for this effect seems to be that the forcing 
scale inputs shear directly into the small scale fields, so the power on small scales 
is the result of cross field structure. But for sufficiently helical turbulence, the 
spectrum changes: the peak at the resistive scale migrates to the forcing scale 
(wavenumber kf) Can we understand the migration of the small scale field 
peak as a function of time? Will this picture survive future numerical testing? 
How do the boundary conditions affect the results? 

To explain the change in the small-scale spectrum from the non-helical to 
helical case, it is possible that what works for the large scale field may also help 
understand what happens for the small scale field. The kinetic helicity input at fc/ 
also cascades to higher wavenumbers, and so there is a source of helicity at these 
wavenumbers. Perhaps the change in shape of the small scale spectrum might be 
modeled by a self-similar set of nested "mean-field" dynamos. The principle is 
that for each small scale wavenumber kg > kf there is a range oi kf < ki < kg for 
which of inverse cascade field growth driven by the helical turbulence at kg can 
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overcome the forward cascade of magnetic energy to fc > fc; . An inverse cascade 
modeled in this way might account for the overall spectral shape change, but 
this is presently just speculation. 

3.2 Large Scale Periodic Box Dynamo 

The critical value of fractional helicity which determines the migration of the 
small scale peak is exactly the same as that for which the kinematic dynamo 
has a positive growth rate, and so a large scale field grows at fci < fc/ in concur- 
rence with the migration of the small scale peak to fc/ (see section 1 and 4). The 
rate of growth, the saturation level, and the dependence on Rm observed in pe- 
riodic box simulations of fully non-linear dynamo seem to be well modeled by 
a two-scale dynamical quenching model of Ref. and section 4. The dynamo 
a is the difference between the kinetic and current helicities associated with the 
small scale field, and the growth of the large scale field is associated with the 
growth of the large scale magnetic helicity. Because total magnetic helicity is 
conserved up to a resistive term, the growth of large scale magnetic helicity also 
means growth of the opposite sign of small scale magnetic helicity, and thus small 
scale current helicity. This eventually suppresses the growth of the large scale 
magnetic field by reducing a. The time dependent process is non-linear. There 
are some unresolved issues with this theory however: 1) the theory works well 
for maximally helical forcing. How does it generalize to arbitrary helical forc- 
ing? 2) How is the correction to a to be properly derived in configuration space 
(see section 5.)? The success of the theory is based on the formula of a from 
which leaves open the ordering ambiguity discussed in section 5. 3) What is 
the appropriate theory for [3 quenching? The saturation value of the mean field 
in the dynamo depends somewhat on the choice of /3 but the ratio of large 
to small scale fields does not. Dynamos with shear and their cycle periods are 
sensitive to the form of (3 4) Along these lines, how does the dynamical 
quenching theory based on magnetic helicity conservation apply to the a — Q 
dynamo? This is work in progress p6| . 5) Do dynamos in astrophysics really 
need to be fast? This is another important reason to study /3 quenching. But 
the question should be revisited. Perhaps the kinematic phase (see (p^) can last 
long enough for substantial amplification, even for asymptotically slow dynamos 
when no cycle period is required. But one must also consider the observed cycle 
periods like that of the sun which seem to be fast. Perhaps highly anisotropic 
turbulence and small cross field structures that can dissipate quickly might even 
allow fast cycle periods with through redefined Rm ■ See also . 

3.3 Coronal Activity and Open Boundary Dynamos (OBD) 

In addition to these considerations, the role of boundary terms needs to be 
considered. In a real system there is gravity, rotation, shear, stratification and 
buoyancy. It may not be enough to just invoke arbitrarily open boundaries to 
test the OBD but the physics at the boundary itself needs to be studied (e.g. 
[^). Winds might be essential. Do mean field dynamos act in symbiosis with 
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jets and winds? We know that the sun is shedding magnetic hehcity. Does this 
enable the mean field dynamo to be fast for all times? In astrophysical rotators, 
the a efi^ect should depend on height. How does this enter in the non-linear 
theory? 

I showed above that the estimated energy deposition rates to maintain fast 
mean field dynamo action are consistent with the coronal + wind power from 
the Sun, Galaxy and Seyfert I AGN all of these sources are natural sites for 
a — i? type dynamos. The helical properties also seem to agree well in the solar 
case where they can be observed |3^. The steady fiow of magnetic energy into 
coronae thus provides an interesting connection between mean field dynamos and 
coronal dissipation in a range of sources. A reasonably steady (over time scales 
long compared to turbulent turnover time scales), active corona with multi- 
scale helical structures, may provide a self-consistency check for fast dynamo 
production of magnetic field. 

For the case, figure 9 shows the effect of including an additional loss term 
in (po|), proportional to /12, on the growth of hi. The effect is to suppress the 
backreaction by taking away excess small scale helicity and allowing the large 
scale field to grow stronger. Understanding appropriate form of loss term for real 
a — fi dynamos needs more study. 



3.4 New Diagram of MFD Operation is Needed 

Related to the boundary issue is the generation of magnetic flux inside of a ro- 
tator. For the Galaxy, diffusion of B at the boundary is required to maintain a 
quadrapole field structure with a net flux inside the disk [Q. The total flux is 
conserved, but to obtain a net flux inside the rotator where desired, the return 
flux must be shed from this region. The surface diffusion of total magnetic field 
for the Galaxy may be difficult and remains an open question. (One alter- 
native possibility is that the mean field diffuses radially more efficiently at the 
top of the disk than at the mid-plane.) Note however, that it is B which needs 
to diffuse, not necessarily the total field or the matter. For the Sun, the solar 
cycle also requires diffusion of mean field through the boundary at the surface 
and/or mid-plane. 

Since both flux generation inside the rotator, and helicity shedding can both 
appeal to boundary terms it is possible that the two are related in the simplest 
a ~ f2 dynamo. Note that it was the growth of the small scale current helicity 
(j-b) that leads to the suppression of a. Thus one can ask directly if the textbook 
a — f2 dynamo has any way of getting rid of this quantity and at the same time 
generating magnetic flux. 

In fact there there is something missing when one draws a standard picture 
of dynamo action either for or for a — f2. The latter is shown in figure 10, 
but the implication holds also shown for the case. The issue is that the small 
scale loop that arises from the kinetic part of the a effect actually induces the 
same sign of the current helicity of the small scale field as that of the large scale 
field that it generates. But I have discussed how the growth of the large scale 
field should be accompanied by the opposite sign of the small scale helicities. It 
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is intriguing because in the kinematic regime of the o? dynamo, the small scale 
field does grow first at the smallest scales, and it seems to take nearly until tkin 
for the peak of the current helicity of opposite sign to the large scale to move 
fully up to the forcing scale (Brandenburg and Maron, personal communication.) 
In this sense, figure 10 seems limited to the kinematic regime: it would be nice to 
see a more accurate picture that actually shows the generation of the small scale 
current helicity of opposite sign to that of the large scale and how it migrates to 
the forcing scale graphically. 

hi 
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250 500 750 1000 1250 1500 1750 2000 

Fig. 9. The effect of a adding a loss term proportional to /12 in (^) on the growth of h\ 
for three different values of loss. The loss term is added by introducing multiplicative 
factors on the last term of (jH^). From top to bottom, these factors are 10, 5 and l(=no 
loss), respectively. Here Rm ~ 200 and k2/k\ — 5. 

Finally, note that it is in part the potential for flux generation that distin- 
guishes a mean field dynamo from the kind of dynamo which generates large 
scale fields as a result of the magneto-rotational instability (MRI). There, the 
field induces turbulence, which then generates a random component of the field, 
which is subsequently sheared by the differential rotation. Magnetic energy grows 
exponentially, even on the largest scale allowed, but there is no real flux gener- 
ation here and no need for input helicity. Note however that for real accretion 
discs, helicity is undoubtedly present since all the ingredients are there, density 
gradient, and turbulence. Thus it is important to understand even in accretion 
disks, what is the role of helical dynamos. This has not yet been done exhaus- 
tively. From this point of view, the mean field dynamo formalism should still 
apply, with a dynamical theory for quenching, and the MRI may play the role 
as the source of turbulence. Nevertheless more work is needed. 
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Fig. 10. The kinematic picture and its limitations. This is assumed to be a northern 
hemisphere so that kinetic helicity is negative, (a) The large scale loop results from 
shear and the small loop results from the a effect. Note that the current helicity of the 
smaller loop has the same sign as that of the large scale loop. So the question arises: 
in a revised picture how can one incorporate the growth of magnetic and/or current 
helicity of the opposite sign that accompanies growth of large scale helicity? (b) This 
just shows the effect of multiple small scale loops as they sum together, (c) the process 
of reinforcement of the toroidal field is shown. 

3.5 What is the Role of Magnetic Reconnection? 

It is generally perceived that reconnection is important for large scale dynamos 
but the precise way in which this is the case is subtle. In general, the processes 
of reconnection can serve two roles: it changes the topology of the field, and also 
dissipates magnetic energy. In the mean field dynamo formalism, the large scale 
field is degenerate with respect to small scale topology. By definition, taking the 
mean means smoothing out over the small scales. Thus a series of disconnected 
loops could have the same mean field as an undulating topologically connected 
field. In the mean field formalism therefore, the role of reconnection is not ex- 
phcit. 

Neither the mean field nor the fluctuating field are the actual field. If one 
really does want the mean field dynamo to result in a topologically connected 
actual field, then reconnection would be important for the topology. But it is 
important to asses the particular application for when this is necessary: the large 
scale field in the Galaxy is measured mainly by Faraday rotation which provides 
little information about the actual field topology. The role of reconnection may 
only be one of ensuring that there is a turbulent cascade which ensures that heat 
rather than the magnetic field is sink of the kinetic energy. 
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Another unanswered question relating to dynamos and rcconncction is: does 
reconnection play in the conservation of magnetic helicity and in the migration 
of the peak of the small scale magnetic energy from the resistive scale to the 
forcing scale in helically forced turbulence? 
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